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ABSTRACT 

Open  loop  frequency  reeponee  measuremente  were  made  of  the 
element!  comprising  the  Foster«Wheeler  DLiG-6^  Class  main  boiler 
and  its  several  control  loops.  Bede  diagrams  were  curve  fitted  and 
the  approximate  transfer  functions  obtained  in  mathematical  form. 
Results  were  analyzed  both  for  agreement  of  open-loop  transfer  functions 
with  continuity  of  energy  and  mass  theory,  and  for  closed  loop  con¬ 
formance  with  stability  criteria.  Data  indicated  that  steam  boiler 
systems  do  not  lend  themselves  readily  to  high  performance  control, 
principally  due  to  large  phase  lags  inherent  in  the  process  at  low 
frequencies.  The  chief  offenders  in  this  regard  are  the  forced  draft 
blowers,  the  heat  sink,  amd  the  additional  integrations  in  the  controllers 


themselves. 
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SUMMARY  PAGE 

The  Plrobleni : 

With  the  advent  of  automatic  combustion  and  feedwater  control 
systems  and  an  ever  increasing  number  of  applications,  to  naval  steam 
generators,  it  has  become  recognized  that  system  designs  might  be 
improved  of  the  mathematical  nature  if  the  systems  were  better  un¬ 
derstood.  As  a  result,  the  Laboratory  has  engaged  in  the  experi¬ 
mental  determination  of  the  dynamic  parameters  of  several  typical 
naval  systems.  This  report  is  concerned  with  the  steam  generator 
system  for  DLG-6  Class, 

Findings: 

The  open  loop  transfer  functions  comprising  the  elements  of  the 
DLG-6  boiler,  blowers,  and  controls,  were  determined  from  measure¬ 
ment  of  their  frequency  response  characteristics. 

Recommendations; 

It  is  recommended  that  the  results  contained  herein  be  made 
available  to  the  systems  industry  in  the  interests  of  enhancing  the 
state  of  the  art  of  control  system  design  for  naval  applications. 
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ADMINISTRATIVE  INFORMATION 
Experimental  determination  of  the  open-loop  dynamic  response 
characteristics  of  the  CVA63,  DLG-6,  and  DLG-9  main  boilers  and  . 
their  control  systems  was  authorized  by  Bureau  of  Ships  letter  NP/4 
Ser  651A  -  993  of  31  July  1961.  The  research  and  development 
pi  oject  number  assigned  was  F013-06-06,  Task  Number  4182,  and 
Laboratory  Project  Number  B-502.  This  is  the  first  of  a  series  of 
final  research  project  reports. 
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REPORT  OF  INVESTIGATION 

INTRODUCTION 

Experience  at  the  Laboratory  and  by  forces  afloat  in  the  operation 
of  automatic  combustion  and  feedwater  control  systems  has  indicated 
that  performance  of  these  systems  is  not  always  in  accordance  with 
specification  requirements  or  designer's  intention.  Such  deficiencies 
generally  result  from  incompatibility  between  the  control  system  and 
the  process  under  control.  This  problem  is  not  easily  overcome  in 
the  design  stages  owing  to  the  lack  of  suitable  information  about  the 
dynamics  of  the  process.  In  an  effort  to  enlarge  the  store  of  in¬ 
formation  relative  to  dynamic  properties  of  systems  in  naval  engi¬ 
neering  plants  the  Bureau  of  Ships  has  authorized  the  Laboratory  to 
conduct  experimental  studies  to  obtain  such  information.  This  report 
covers  the  investigation  and  results  obtained  in  one  such  study,  the 
OLG-6  main  boiler,  blowers,  and  automatic  controls. 

OBJECT 

The  object  of  this  project  was  to  measure  the  frequency  or  transient 
response  characteristics  in  open  loop  of  the  elements  comprising  the 
DLG-6  Class  main  boiler  and  its  accessories,  and  to  express  these 
characteristics  mathematically  in  the  form  of  open-loop  transfer  functions. 


1 


NBTL  PROJECT  B-SOZ-I 


The  frequency  response  technique  was  to  be  used  for  naeasurement 
wherever  applicable. 

PROCEDURE 

The  Laboratory's  DLG-6  boiler  Installation  (NBTL  Project 
B-27P)  was  equipped  with  instrumentation  to  provide  electrical 
signals  proportional  to  the  following  variables: 

1.  Combustion  and  Feedwater  Controls 

a.  Steam  pressure  controller  output 

b.  Air  flow  controller  output 

c.  Selective  relay  output 

d.  Characterizing  relay  output 

e.  Feedwater  controller  output 

f.  Air  flow  transmitter  output 

g.  Steam  pressure  error 

h.  Air  flow  error 

2.  Boiler 

a.  Steam  flow 

b.  Feedwater  flow 

c.  Return  fuel  oil  pressure 

d.  Water  level 
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e.  Final  ateam  temperature 

f.  Superheater  outlet  presaure 

3.  Forced  Draft  Blowera 

a.  Noazle  valve  lifting  gear  poaition 

b.  Forced  draft  blower  apeed 

An  electronic  frequency  reaponae  analyzer,  Boonahaft  U  Fucha, 

Inc. ,  Type  71  lA,  was  utilized  to  measure  phase  and  amplitude  re¬ 
lationships  between  sinusoidal  variations  of  pairs  of  variables  as  listed 
above.  This  device  was  equipped  with  an  oscillator  (adjustable  from  .01 
to  1000  cycles  per  second)  for  excitation  of  the  system  through  an 
electro-pneumatic  converter.  The  analyzer  is  desigpied  to  compare 
only  the  fundamentals  at  the  desired  frequencies,  rejecting  40  db  of 
all  other  components;  this  feature  permits  a  true  frequency  analysis 
despite  the  introduction  of  non-linear  distortions  elsewhere  in  the 
system.  In  cases  where  a  frequency  response  method  was  not  applicable, 
the  open-loop  transient  response  characteristics  were  obtained  by 
recording  responses  to  step  input  disturbances  using  a  two-channel 
recording  oscillograph.  In  cases  where  the  range  of  interest  of  fre¬ 
quencies  fell  below  the  capabilities  of  the  analyzer,  a  low  frequency 
oscillator  and  the  recording  oscillograph  were  employed. 

Data  were  obtained  by  the  above  techniques  and  plotted  in  the  form 
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of  Bode  Diagrams.  These  diagrams  were  subsequently  analyzed  to 
determine  the  mathematical  open-loop  transfer  functions. 

It  will  be  noted  that  the  procedure  outlined  above  produces  the 
perturbed  form  of  the  transfer  function  linearized  in  a  narrow  range 
about  the  mean  operating  level.  This  is  not  descriptive  of  the  system 
non-linearities,  but  does  provide  a  reasonable  basis  for  system  design 
or  analysis.  The  nature  and  extent  of  these  non-linearities  are  to  be 
investigated  following  major  revisions  to  the  electronic  equipment 
employed  for  this  type  of  work. 

RESULTS 

The  results  of  this  project  are  shown  graphically  in  Figures  1 
through  19.  Taken  one  element  at  a  time,  the  closed  loop  system 
may  be  described  by  the  following  constituent  parts: 

1.  Automatic  Combustion  Controls 

a.  Steam  pressure  controller:  This  unit,  a  Moore  Products 
Co.  Model  5  >  Proportional  plus  Reset  stacked  diaphragm  type  controller, 
was  found  by  its  frequency  and  transient  response  characteristics  to  be 
in  unsatisfactory  working  condition.  The  Bode  diagram  for  this  d'^vice, 
shown  in  Figure  1,  does  not  conform  to  the  usual  proportional  plus  re¬ 
set  characteristics.  The  ideal  response  characteristic  curves  for  this 
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initrument  are  indicated  by  the  solid  lines  of  Figure  1.  It  is 
apparent  from  the  experimental  data  that  the  reset  mode  is  not  present 
at  low  frequencies.  This  of  course,  is  quite  the  opposite  effect  of 
that  anticipated,  as  indicated  by  the  solid-line  curves.  The  idekl 
curves  are  located  with  ruspect  to  the  gain  and  frequency  scales 
according  to  the  values  of  gain  and  reset  time  previously  determined 
during  control  system  tests.  In  an  effort  to  determine  the  cause  of 
the  unusual  Bode  diagram,  the  controller  was  dismantled,  inspected, 
cleaned,  and  reassembled.  Retests  produced  essentially  identical 
results.  A  transient  response  test  in  which  a  step  change  was  made 
in  the  error  signal  produced  a  response  curve  resembling  proportional 
plus  exponential.  In  service,  the  unit  provided  control,  but  with  some 
proportional  offset  and  insufficient  recovery.  As  an  estimate  of  the 
performance  of  a  properly  adjusted  controller,  the  transfer  function 
shown  in  Figure  1  was  obtained  from  the  idealized  Bode  curves  of 
Figure  1. 

b.  Air  flow  controller:  Because  of  the  relatively  low 
natural  frequency  of  this  device,  a  Moore  Products  Co.  Model  68VT17 
Proportional  plus  Reset  Controller,  it  was  initially  decided  to  measure 
the  open  loop  response  to  a  step  input  to  determine  the  proportional 
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gain  and  reset  rate.  The  step  changes  in  the  input  signal  were  made 
in  the  upward  and  downward  directions  with  the  following  results: 

Step  Increase:  Gain  =  1. 23;  Reset  rate  »  .6lrep/min 

Step  Decrease:  Gain  >  1.37;  Reset  rate  ■  .52  rep/min 

Dead  time  ■  .002  sec 

The  average  transfer  function  can  therefore  be  expressed  as: 

G(t)  =|l.3  + 

It  is  noted  that,  in  contrast  to  the  steam  pressure  controller, 
this  device  exhibited  the  customary  proportional  plus  integral  transient 
response.  The  magnitude  of  the  step  changes  was  0.7  psi,  starting 
with  an  initial  9  psi  input  and  output. 

The  needle  valve  dial  on  this  instrument  indicated  a  nominal 
setting  of  .23  repeats  per  minute  reset  rate;  the  gain  is  of  course  non- 
adjustable  and  is  designed  for  unity.  Indications  are  that  both  of  these 
values  are  substantially  in  error. 

c.  Fuel  oil  burner  characterizing  relay;  This  device, 
General  Regulator  Co.  Model  1771,  is  designed  to  convert  a  controller 
output  representing  demand  for  fuel  oil  flow  per  burner  to  a  control 
valve  loading  pressure  representing  demand  for  fuel  oil  pressure. 
Characterizing  is  done  by  means  of  a  cam  arrangement  programmed  to 
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provide  the  proper  relationship  between  these  functions  as  dictated 
by  the  pressure/flow  characteristics  of  the  particular  sprayer  plates 
selected  for  the  fuel  oil  burners.  For  OLG-6  class,  this  arbitrary 
function  program  of  the  characterising  relay  is  illustrated  in  Figure  2. 

The  frequency  response  of  the  actuator  was  measured  and  the 
resulting  Bode  Diagram  is  shown  in  Figure  3.  The  solid  line  indicates 
experimental  values,  while  the  dotted  line  indicates  the  approximation 
obtained  by  plotting  the  transfer  function  shown.  The  gain  scale  was 
located  such  that  zero  db  represents  an  amplitude  ratio  of  2.08 
psi/psi,  the  zero  frequency  gain  measured  at  the  boiler  cruising 
condition.  The  gain  is  of  course  variable  with  operating  level,  and 
can  be  obtained  from  the  slope  of  the  curve  shown  in  Figure  2  for  any 
average  value  of  input. 

The  transfer  function  is  of  the  usual  quadratic  form  associated 
with  the  second  order  response  of  servo-mechanisms  of  this  type. 

The  transfer  function  is  factorable  into  two  series-connected  first 
order  linear  lags  having  time  constants  of  0.41  and  1.12  seconds. 

The  existence  of  real  roots  indicates  considerable  overdamping.  The 
Bode  diagram  also  approximates  a  pure  second  order  system  with 
overdamping  and  a  natural  frequency  of  0.24  cycles  per  second. 
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d.  Force 'balance  fuel  oil  control  valve:  Thia  unit  responds 
to  the  output  of  the  characterizing  relay  to  produce  corresponding 
back  pressure  in  the  return  fuel  oil  burner  line.  It  is  hydraulically 
powered  by  the  return  fuel  oil  pressure.  Results  of  frequency  response 
testa,  again  at  the  cruising  condition,  are  shown  in  Bode  diagram  form 
in  Figure  4.  The  transfer  function  is  again  a  factorable  lag,  with  time 
constants  of  each  component  0.  355  and  0.455  seconds.  The  continued 
phase  shift  at  higher  frequencies  indicates  that  a  small  amount  of 
backlash  exists;  the  dead  time  associated  with  this  backlash  amounts 

to  0.023  seconds.  The  zero  frequency  gain  is  15  db,  or  5.6  psi  oil 
pressure  per  unit  psi  pneumatic  loading  pressure. 

e.  Forced  draft  blower  actuator:  The  Model  1801  General 
Regulator  Co.  pneumatic  power  actuator  positions  the  throttle  linkage 
of  the  Hardie-Tynes  main  forced  draft  blowers  amd  the  combustion 
control  dampers.  Frequency  response  tests  were  conducted  by 
measuring  the  angular  displacement  of  the  drive  crank  with  respect 
to  sinusoidal  variations  in  the  input  position  demand  signal  about  a 
fixed  reference  corresponding  to  the  normal  position  with  two  blowers 
in  operation  at  the  cruising  condition.  ,  The  Bode  diagram,  Figure  5, 
is  a  typical  second  order  linear  underdamped  response.  The  transfer 
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function  hni  the  dimeneions  "degreee  rotation  per  unit  pei  pneumatic 
input  signal",  a  natural  frequency  at  .89  cycles  per  second,  and  a 
damping  factor  of  .  23.  This  transfer  function  is  typical  of  positioning 
servo-mechanisms  of  this  type. 

f.  Air  flow  transmitter:  This  device,  General  Regulator 
Corp.  Model  30244,  senses  the  differential  air  pressure  across  the 
burner  registers,  extracts  the  square  root  of  this  differential  pressure 
and  transmits  a  proportional  pneumatic  output  signal.  Frequency 
response  tests  were  conducted  by  discharging  the  output  of  an  electro¬ 
pneumatic  transducer  through  a  small  orifice,  and  metering  the  drop 
across  the  orifice  with  the  transmitter.  Sinusoidal  excitation  of  the 
e-p  transducer  thus  produces  a  sinusoidal  transmitter  output.  The 
orifice  differential  pressure  was  related  to  the  corresponding 
differential  across  the  registers  with  all  burners  in  use  for  known 
values  of  steady  state  air  flow,  thus  establishing  the  zero  frequency 
gain  of  the  transmitter.  Data  shown  in  the  Bode  diagram.  Figure  6, 
illustrate  the  typical  quadratic  response  of  force -balance  type  in¬ 
struments. 

Curve  fitting  of  the  Bode  diagram  results  in  very  close  agreement 
with  the  second  order  transfer  function  with  .  50  damping  factor.  The 
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zero  frequency  gain  is  measured  in  cubic  feet  of  air  per  minute  to  the 
boiler  furnace,  referred  to  standard  conditions.  This  gain  was  found 
to  be  2.41  X  10**^  psi  output  per  acfm  of  air  flow  to  the  boiler.  The 
natural  frequency  of  the  transnnitter  was  located  at  12. 6  cycles  per 
second. 

g.  Copper  Tubing:  Measurement  was  made  of  the  fre¬ 
quency  response  of  the  1/4"  O.  D.  copper  tubing  connecting  the  ratio 
relay  to  one  forced  draft  blower  actuator,  a  distance  of  approximately 
66  feet,  with  negligible  end  displacement  or  capacitance.  The  Bode 
diagram.  Figure  7,  reveals  the  anticipated  fourth  order  system  with 
natural  frequencies  of  4.0  and  6.  3  cycles  per  second.  It  is  to  be 
noted  that  the  tubing  contributes  a  gain  of  nearly  3  db  at  frequencies 
in  the  range  of  two  to  three  cycles.  This  may  contribute  to  the  high 
frequency  jitter  often  observed  in  blower  actuators  and  air  flow  feed¬ 
back  signals  at  the  control  board.  The  characteristics  of  otlier  long 
runs  of  tubing  can  be  inferred  from  the  experimental  results  for  this 
particular  section. 

2.  Feedwater  Regulator 

a.  Feedwater  control  valve:  The  Copes -Vulcan  Model 
CV-D  5"  Main  Feedwater  Regulating  Valve  was  tested  by  exciting 
the  input  to  the  valve  positioner  and  measuring  the  response  of  feed- 


10 


NBTL  PROJECT  B-SOZ-I 


water  flow  rate  while  maintaining  an  average  differential  preesare 
acroea  the  valve  of  75  pei.  Reeulta,  shown  in  Figure  8,  were  very 
closely  approximated  by  a  third  order  transfer  function  with  back¬ 
lash.  This  transfer  function  consists  of  a  first  and  a  second  order 
component  with  natural  frequencies  at .  20  and  .158  cycles  per 
second  respsctivsly.  The  first  order  lag  is  believed  to  be  that  of 
the  valve  and  diaphragm  operatori  with  the  time  constant  associated 
with  the  capacitance  of  the  diaphragm  chamber:  the  second  order 
lag  is  believed  to  be  descriptive  of  the  Moore-Products  Co.  valve 
positioner  used  with  the  regulating  valve.  It  is  noted  that  considerable 
under  damping  is  present,  with  the  result  a  six  db  gain  at  .  15  cycles 
per  minute.  The  third  component,  .45  seconds  of  dead  time,  is  attribu¬ 
table  to  stick-slip  friction  of  the  valve  stem,  probably  in  the  stuffing 
box.  Some  of  this  dead  time  might  be  eliminated  by  installation  of 
fresh  packing  and  proper  adjustment  of  the  packing  gland. 

3.  Hardie-Tynes  Main  Forced  Draft  Blowers. 

a.  Blower  speed/air  flow  demand:  Variation  in  forced  draft 
blower  speed  in  response  to  excitation  of  the  actuator  input  signal  are 
illustrated  in  Bode  diagram  form  in  Figure  9.  Response  was  measured 
by  perturbation  of  the  normal  input  signal  when  operating  with  two 
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blowers  at  the  cruieing  condition;  combustion  control  dampers  were 
open  under  this  condition.  This  repponse  is  considered  to  be  grossly 
non-linear  due  to  characteristics  of  the  steam  turbine  noszle  valves 
and  the  non-linear  speed-power  relationship.  The  perturbed  form 
of  the  transfer  function  reveals  two  linear  lags,  the  first  associated 
with  the  Inertial  properties  of  the  rotating  system,  and  the  second 
with  the  compressibility  effects  of  the  air  at  the  blower  discharge. 

The  zero  frequency  gain  varies  with  the  operating  level,  and  can  be 
obtained  from  the  steady  state  relationship  between  speed  and  actuator 
input  shown  in  Figure  10.  The  dynamic  non-linearity  can  be  expressed 
in  terms  of  blower  speed;  this  relationship  will  be  further  explored 
in  a  subsequent  section  of  the  report. 

b.  Air  flow /blower  speed  :  The  air  flow  is  directly  pro¬ 
portional  to  the  blower  speed  with  all  burner  registers  open,  in 
accordance  with  steady  state  fan  laws.  Considering  the  dynamic 
compressibility  effects  to  be  included  in  the  speed  transfer  function,  . 
the  air  flow/speed  function  can  be  considered  to  be  pure  gain.  A  plot 
of  blower  capacity  vs.  blower  speed  produces  a  linear  function  with 
a  slope  of  4.95  scfm  per  unit  RPM,  for  one  blower.  This  slope  is 
equivalent  to  the  gain  required  to  express  the  transfer  function  obtained 
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in  2(a)  above  in  terms  of  air  flow  per  unit  air  flow  demand  signal. 

4,  Main  Boiler  Energy  and  Mass  Balance. 

The  energy  balance  across  the  boiler  under  dynamic 
conditions  determines  the  magnitude  of  energy  stored  within  the 
boiler,  which  in  turn  fixes  the  output  conditions,  pressure  and  tern • 
perature.  The  open>loop  transfer  function  describing  the  properties 
of  this  heat  balance  may  be  listed  as  follows:  steam  pressure/heat 
input,  steam  pres  sure /heat  output,  steam  temperature /heat  input, 
steam  temperature /heat  output,  water  level/heat  input,  water  level/ 
heat  output,  water  level/water  flow,  and  steam  pressure/water  flow. 

Procedure  for  obtaining  these  transfer  functions  included 
both  transient  and  frequency  response  methods.  In  both  cases  in* 
teractions  were  avoided  insofar  as  was  possible  by  placing  all  para¬ 
meters  in  open-loop,  and  measuring  responses  between  output  and 
input  pairs.  Detailed  results  were  as  follows: 

a.  Steam  pressure/heat  input:  The  response  of  super¬ 
heater  outlet  pressure  to  disturbances  in  heat  input  was  measured 
by  introducing  a  sinusoidal  oscillation  into  the  fuel  rate  through  the 
burners  while  maintaining  constant  steam,  feed,  and  air  flow  rates. 
The  system  was  otherwise  permitted  to  remain  in  closed-loop  and  the 
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excitation  introduced  at  the  set  point.  Although  operation  was  in  closed- 
loop  only  the  open-loop  response  was  measured.  The  Bode  diagram  is 
shown  in  Figure  11.  Curve  fitting  yielded  a  very  good  agreement  with 
the  gain  curve,  but  only  a  fair  approximation  of  the  phase  curve.  The 
approximate  transfer  function  was  expressed  as  two  series-connected 
linear  lags  having  widely  separated  natural  frequencies.  In  general, 
the  high  frequency  component  can  be  neglected  with  no  appreciable 
effect  on  the  mathematical  description  of  the  system. 

The  inability  to  fit  the  phase  curve  with  accuracy  implies 
that  dynamic  non-linearities  exist.  Two  principal  non-linearities  are 
the  heat  absorbtion  characteristics  of  the  generating  section,  and  the 
thermodynamic  properties  of  the  steam -water  mixture  in  the  boiler. 

An  additional  static  non-linearity  results  from  the  fuel  oil  burner 
pressure -flow  characteristics.  This  non-linearity  can  be  eliminated 
by  expressing  the  input  in  heat  rate  or  fuel  rate  terms,  as  indicated  in 
Figure  11. 

b.  Steam  pressure/heat  output:  The  response  of  super¬ 
heater  outlet  pressure  to  changes  in  steam  flow  was  measured  by  the 
transient  response  method.  This  procedure  was  necessitated  by  the 
fact  that  the  Laboratory's  main  steam  throttle  valve  is  driven  by  a  large 
induction  motor.  No  means  were  available  for  producing  oscillatory  steam 
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flow  variations,  but  the  high  stroke  speed  available  with  the  motor  drive 
permitted  a  near  step  change  inwalve  position.  A  step  change  of  minus 
two  percent  of  load  was  made  starting  from  the  cruising  steam  rate. 

The  response  of  superheater  outlet  pressure  was  recorded.  Values 
of  the  incomplete  response  plotted  against  time  on  semi-logarithmic 
paper  plotted  to  nearly  a  straight  line,  indicating  an  exponential,  or 
first  order  linear  lag  type  response.  The  single  time  constant,  the 
value  of  t  at  36.8  percent  incomplete  response,  was  6.4  seconds.  A 
second  component,  producing  linear  decay  with  time,  was  identified  as 
an  integration.  When  plotted  on  rectangular  paper  as  shown  in  Figure  12, 
the  integration  component  exhibited  a  slope  of  .  383  psi/second  for  a  steam 
♦low  step  change  of  4200  Ihs/hr.  From  these  data,  an  approximate 
trajisfer  function  can  be  derived: 


s  =  -  1. 37  X  10 
Ws 


Y  r+'i'T')  ~ -91  X  10”^(-U  ,  _£^ 

V  1  +  6.  4p  y  V  P  /  lbs/] 


SI _ 

/hr 


Where  Pg  =  Superheater  cutlet  pressure,  and  "Wg  -  steam  output  rate. 

c.  Steam  temperature /heat  input:  Superheater  performance 
data  vere  obtained  by  varying  the  fuel  input  rate  and  measuring  the  super¬ 
heater  outlet  temperature  response.  The  steam  temperature  is  a  function 
of  t'^.c  heat  absorbtion.  in  the  superheater,  which  in  turn  results  from 
cond'.- ction,  convection,  and  ladiatior.  in  the  steam  generator.  Analyzed 
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computation  of  ateam  temperature  by  claasical  heat  tranafer  equationa 
reveala  the  preaence  of  groaa  non-linearitiea  in  the  ateady  atate  and 
tranaient  performance.  Aa  ia  to  be  expected,  thia  ia  reflected  in  the 
non-linearity  of  the  frequency  reaponae  characteriatica  obtained 
experimentally.  The  Bode  diagram  of  theae  characteriatica  appeara 
in  Figure  13;  all  attempta  to  curve  fit  linear  equationa  to  theae  curvea 
to  obtain  approximate  tranafer  functiona  were  to  no  avail  becauae  of 
the  unuaual  amplitude -phaae  relationahipa.  Inapection  of  Figure  13 
diacloaea  large  phaae  laga  with  increaaing  frequency  accompanied 
by  relatively  email  reduction  of  gain.  Thia  implies  the  presence 
of  considerable  transport  delay  time.  Transport  delay  time  is  of 
course  a  common  property  of  heat  exchangers;  but  the  high  steam 
velocities  in  the  superheater  tubes  and  piping  lend  small  credence  to 
the  existence  of  substantial  amounts  of  distance -velocity  lag  between 
the  "center  of  heat  absorbtion"  and  the  measuring  element.  On  the 
other  hand,  if  it  is  assumed  that  the  superheater  has  a  very  low  heat 
capacity,  then  variations  in  available  energy  in  the  gas  path  can  be 
expected  to  produce  corresponding  variations  on  the  steam  side.  Changes 
in  phase  can  be  attributed  to  the  heat  transfer  resistance  of  the  steam 
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and  gat  boundary  layer  films.  On  this  basis  it  can  be  argued  that  the 

drop  in  gain  with  increasing  frequency  is  attributable  to  the  first  order 

lag  associated  with  superheater  tube  metal  heat  capacity,  and  the 

additional  phase  shift  to  several  distance -velocity  lag  components  at 

relatively  high,  but  different  frequencies. 

d.  Steam  temperature /heat  output:  The  main  steam  throttle 

valve  was  not  capable  of  sinusoidal  oscillation:  furthermore,  because 

of  the  gross  non-linearity  of  the  superheater  response,  as  discussed 

above,  it  was  recognized  that  results  of  a  transient  response  test 

would  be  essentially  meaningless,  or  at  best  a  very  rough  approximation. 

Accordingly,  an  analytical  determination  of  the  response  was  considered 

to  be  more  truly  descriptive  of  performance  than  that  to  be  obtained 

by  any  experimental  method  available.  The  steady  state  energy  balance 

across  the  superheater  may  be  written  thus: 

Wa  (hg  -hg)  -  K  Wo  (HHV) 

Where:  Wg  •  Steam  flow,  lbs /hr 

hg  ai  Enthalpy  of  superheated  steam,  Btu/lb 

hg  -  Enthalpy  of  saturated  steam,  Btu/lb 

K  «  Fraction  of  useful  heat  absorbtion  from 
combustion  credited  to  the  superheater 
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Wq  b  Fuel  ratej  lbs /hr 

HHV  <■  Higher  heating  valire  of  fuel,  Btu/lb 

e^  B  Boiler  efficiency 

Subatitution  of  numerical  valuea  appropriate  to  the  cruiaing 
condition  givea: 

56.000  (1462  -  1183)  «  K(4130)  (18,  500)(.  878) 

Solving  for.l^, 

(56.000)  (279)  _ 

(4130)  (18,  500)  (.878)  "" 

Rewriting,  Wg  (279)  =  Wq(3740) 


or 

J£o  =  .075 
Ws 

It  can  now  be  assumed  that  the  Bode  diagram  of  steam  temperatur^ 
steam  flow  is  identical  to  that  of  Figure  13,  with  the  gain  scale  shifted 
by  a  factor  of  .075,  or  -22.5  db. 

e.  Water  ievel/heat  input;  The  frequency  response  of  boiler 
drum  water  level  with  respect  to  disturbances  in  heat  input  was 
determined  by  oscillating  the  fuel  oil  control  valve.  A  Bode  plot  of  the 
resulting  data  is  shown  in  Figure  14.  Analysis  of  the  Bode  diagram 
indicated  that  the  transfer  function  could  be  closely  approximated  by 
two  series -connected  first-order  linear  lags,  with  natural  frequencies 
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at  .063  and  .279  cycles  per  second.  It  is  likely  that  these  two  com¬ 
ponents  represent  the  volume  displacement  changes  with  respect  to 
heat  abaorbtion  changes,  and  mass  displacement  changes  with 
respect  to  changes  in  differential  head  (hydraulic  lag),  respectively. 

f.  Water  level/heat  output:  Response  of  water  level  to 
changes  in  steam  flow  were  measured  by  the  transient  response  method, 
again  due  to  inability  to  oscillate  the  steam  valve.  Step  chamges  in 
steam  flow  at  the  cruising  condition,  and  the  resulting  variations  in 
water  level  (holding  water  flow  and  oil  flow  constant)  Were  recorded. 

The  transient  response  is  shown  in  Figure  15.  From  these  transient 
response  data  the  transfer  function  was  recognized  to  consist  of  two 
parallel  components:  a  linear  lag,  and  an  integration.  The  first  of  these 
is  attributable  to  energy  considerations  resulting  from  chatnges  in  heat 
absorbtion  in  the  generating  tubes.  This  portion  of  the  transient  response 
was  further  investigated  by  plotting  incomplete  response  against  time 
on  semi -logarithmic  paper,  and  obtaining  the  time  constant  at  37.8 
percent  incomplete  response.  The  average  value  for  the  upward  and 
downward  step  changes  was  found  to  be  4.  3  seconds.  The  second 
component,  the  integration,  results  from  continuity  of  mass  consideration. 
A  weight  rate  balance  across  the  boiler  indicates  that  the  integration 
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result!  from  a  difference  between  steam  flow  and  water  flow,  and  that 
the  integration  rate  is  proportional  to  this  difference  and  dependent  upon 
the  boiler  drum  geometry.  This  relationship  is  further  discussed  in 
the  following  section. 

g.  Water  level/water  flow :  Sinusoidal  variation  in  feed- 
water  flow  about  the  normal  cruising  condition  was  generated  by  osci¬ 
llating  the  main  feed  regulating  valve  through  a  series  of  test  frequencies. 
The  response  of  water  level  in  the  steam  drum  to  these  disturbances 

was  measured.  The  resulting  frequency  response  Bode  diagram  appears 
in  Figure  16.  Analysis  of  the  diagram  revealed  that  the  transfer  function 
consisted  of  two  components:  a  linear  lag  auid  an  integration,  series - 
connected.  Despite  the  fact  that  the  range  of  test  frequencies  was 
considerably  higher  than  the  natural  frequency  of  the  low  frequency 
component,  it  was  readily  identified  as  a  true  integration  from  the  con¬ 
tinuity  of  mass  equation,  in  accordance  with  the  same  argument  ad¬ 
vanced  in  (f)  above. 

h.  Steam  pressure/water  flow:  Efforts  to  obtain  the  response 
of  steam  pressure  to  changes  in  water  flow  revealed  that  the  frequencies 
of  interest  were  well  below  the  capabilities  of  the  frequency  response 
analyzer;  attenuation  of  the  steam  pressure  signal  at  higher  frequencies 
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produced  considerable  scatter  in  the  experimental  data.  The  transfer 
function  can,,  however,  be  deduced  from  energy  consideration.  Consider 
the  energy  balance: 

h  .  k/(KiW„  -  K,W.  .  K,W„  .  K,AW„)dt 

Where:  hf  «  change  in  enthalpy  of  saturated 
liquid,  Btu/lb 

Wq  *  Fuel  rate,  lbs /hr 

Ws  *  Steam  rate,  lbs /hr 

W^.m  Water  rate,  Ibs/hr 

Wm  *  Weight  of  water /steam  mixture  in 
the  boiler,  lbs 

K  •  Constant 

The  first  three  terms  in  the  integral  represent  steady  state  con¬ 
ditions,  and  their  sum  must  be  equal  to  zero.  The  value  of  the  integral 
when  divided  by  the  weight  of  the  contents  of  the  steam  generating 
system  is  a  measure  of  the  unit  energy  of  the  boiler  mixture  and  roughly 
equal  to  the  enthalpy  of  the  saturated  liquid  at  the  operating  pressure. 

For  a  small  perturbation  in  water  flow  the  weight  of  boiler  contents 
can  be  considered  constant.  Therefore  the  equation  can  be  rewritten: 
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In  view  of  the  fact  that  the  saturated  liquid  enthalpy  varies 
linearly  with  pressure  in  the  range  of  operating  pressures  under 
consideration,  the  above  maybe  expressed  as: 


Or  in  transfer  function  form: 


The  same  analysis  may  be  applied  to  the  steam  flow  and  fuel  flow 
terms  of  the  basic  equation.  The  natural  frequency  in  each  case  is  a 
function  of  boiler  geometry,  capacity,  ^circulation  and  heat  abosorbtion 
characteristics,  and  operating  pressure.  Transfer  functions.for  steam 
flow  or  feedwater  flow  are  related  by  the  ratio:  of  K3  to  K^,  which  in 
turn  are  proportional  to  the  heat  absorbtion  in  boiling  and  in  pre¬ 
heating  to  saturation.  At  the  cruising  condition: 


-  .Af  -  215-366 


608 


and 


(fn)  =  .25(f„) 

feedwater  steam 

flow  flow 


From  the  integration  term  of  the  steam  pressure /steam  flow 
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transfer  function  of  Figure  12,  the  transfer  function  being  sought 
(neglecting  low  frequency  effects)  is: 


This  transfer  function  is  analytically  derived,  but  based  on  the 
experimentally  determined  value  of  the  natural  frequency  of  a  related 
function. 

CONCLUSIONS  AND  DISCUSSION 

Conclusions  pertaining  to  results  previously  described  are  in 
corresponding  order; 

1.  Automatic  Combustion  Controls: 

a.  The  steam  pressure  controller  did  not  exhibit  the  an¬ 
ticipated  proportional  plus  reset  characteristics,  thus  explaining 
subsequent  marginal  performance  of  the  control  system  when  operating 
in  closed  loop.  Inspection  of  this  component  failed  to  disclose  the 
reason  for  its  unusual  characteiistic  response.  In  general,  it  has  been 
the  Laboratoiy's  experience  that  the  stack  diaphragm  type  controllers, 
particularly  those  employing  pressure-splitting  gain  circuitry  exhibit 
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troublesome  unreliable  performance.  In  any  case,  no  credence  should 
be  given  to  the  dial  readings  of  calibrated  needle  valves  used  in  in¬ 
struments  of  this  type. 

b.  Air  flow  controller  taansient  response  conformed  to 
anticipated,  in  contrast  to  that  of  the  steam  pressure  controller. 
Non-linearity  of  this  device  was  revealed  in  the  varying  gain  and 
reset  rate  values  obtainedi  under  increasing  and  decreasing  step  input 
conditions.  Again  the  reset  rate  needle  valve  dial  did  not  agree  with 
the  measured  reset  value. 

c.  The  fuel  oil  characterizing  relay  is  essentially  a 
pneumatic  force-balance  servo-mechanism.  The  quadratic  lag  measured 
in  frequency  response  tests  was  considered  typical  of  devices  of  this 
type. 

d.  The  fuel  oil  control  valve  also  exhibited  a  characteristic 
quadratic  lag  with  considerable  overdamping.  It  is  interesting  to  note 
that  this  device  has  a  better  frequency  response  than  the  characterizing 
relay  which  drives  into  it. 

e.  The  response  of  the  forced  draft  blower  pneumatic 
throttle  actuator  was  considered  typical  of  the  characteristics  of  units 
of  this  type. 
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f.  The  response  of  the  air  flow  transmitter  was  again  a 
typical  second  order  function  usually  associated  with  force -balance 
instruments. 

g.  The  characteristic  Bode  diagram  for  copper  tubing  was 
of  a  form  coincidental  with  published  results  of  other  investigations. 

The  analysis  of  the  dynamics  of  transmission  lines  is  extremely  complex 
and  is  generally  considered  too  complicated  for  practical  engineering 
applications.  Assuming  that  the  receiver  volume  remain  small  in 
relation  to  the  volume  of  the  tubing,  use  of  an  inverse  proportion 
relationship  between  tubing  length  and  natural  frequency  for  each  of 

the  quadratic  terms  will  provide  a  useful  approximation.  The  transfer 
function  indicated  in  Figure  7  is  for  a  section  of  dead  ended  tubing,  66  feet 
in  length.  Other  transmission  lines  of  roughly  equal  lengths  are  found  at 
the  output  of  the  air  flow  transmitter  and  steam  pressure  transmitter. 

2.  Feedwater  Regulator 

a.  The  feedwater  control  valve  characteristics  conform 
closely  to  a  second  order  lag  with  additional  first  order  and  backlash 
components  representative  of  diaphragm  chamber  capacity  and  friction. 
The  backlash  dead  time  of  one  half  second  is  believed  to  be  excessive, 
but  fortunately  does  not  appear  to  introduce  excessive  phase  shift  in  the 
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closed  loop.  This  backlash  mij^t  be  relieved  by  packing,  lubrication, 
or  similar  maintenance. 

3.  Main  Forced  Draft  Blowers. 

a.  The  Bode  diagram  of  Figure  9  illustrates  the  per¬ 
turbed  frequency  response  of  one  blower  operating  at  the  two-blower 
cruising  condition.  Recognizing  that  the  transfer  function  is  non¬ 
linear,  the  experimental  results  may  be  compared  to  analytical  results 
as  follows:  The  angular  acceleration  of  the  blower  rotor  is  directly 
proportional  to  the  net  applied  torque,  and  inversely  proportional  to 
the  polar  moment  of  inertia.  Hence: 


dtp  _1 _  ^ 

dt  *  MK^  ’ 


(Ti  -  T2) 


Where:  <»l  •  Angular  velocity 

M  *  Mass  of  rotating  system 
K  s  Radius  of  gyration 
T  j  •  Torque  applied  by  turbine 
T2  a  Torque  applied  by  fan 

In  terms  of  revolutions  per  minute  and  poler  moment  of  inertia, 
and  integrating: 
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Where:  N  «  Angular  velocity,.  RPM 

Ip  *  Polar  moment  of  Inertia 
K  «  Constant 

p  ■  Differential  operator, 

The  torque  applied  by  the  fan  acts  as  a  brake;  this  torque  varies 
directly  with  the  square  of  angular  velocity.  Therefore: 

i4(Ti-CN“?)  _L_ 

M  Ip  •*  p 

„4  4..K,  (i-  £nAj. 

T,  IV  Tj/p 

This  equation  is  the  transfer  function  of  the  response  of  blower 
speed  to  changes  in  applied  torque .  In  perturbed  form  the  applied 
torque  can  be  considered  proportional  to  steam  flow,  and  the  braking 
torque  proportional  to  rotational  speed. 

2K  XT  1  1 

N._  (CjW,-CjN)- 
P  ^ 


NeCjWg 


C2N 


1  + 


60  C 


2TTI, 


2)«C,Wa 


P 


27 


NBTL  PROJECT  B -502-1 


N  60C.(l4>) 


9.55C,  _  C,/C>  ^ 

•  Ip(p+9.55^^  (1+  .1047(-^)  P) 

The  transfer  function  now  appears  in  first  order  linear  lag  form, 
where: 

N  •  Rotational  speed,  RPM 

Wg  c  Steam  flow  to  turbine 

C^«i  Slope  of  torque/steam  flow  curve 

C2«  Slope  of  torque/speed  curve 

The  time  constant  has  the  value  .  1047  ,  indicating  that  rapid 

dynamic  response  can  be  expected  of  blower  designs  having  minimum 
polar  moments  of  inertia  and  maximum  rate  of  change  of  torque  with 
respect  to  angular  velocity.  The  transfer  function  derived  above  con¬ 
stitutes  the  principal  dynamic  lag  in  the  response  of  the  forced  draft 
blower,  and  is  equivalent  to  the  lower  frequency  lag  components  of  the 
experimentally  determined  transfer  function  shown  in  Figure  9.  The 
higher  frequency  component  probably  arises  from  the  resistance  and 
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capacitance  at  the  blower  discharge.  This  RC  network  will  result  in 
a  phase  lag  between  speed  and  air  flow,  which  will  produce  a  corre¬ 
sponding  lag  between  speed  and  static  pressure,  hence  speed  and 
power,  hence  speed  and  retarding  torque.  Because  of  the  complex 
s/stem  geometry  the  numerical  value  of  this  capacitance  is  difficult 
to  compute  with  any  degree  of  accuracy.  In  general  it  appears  that 
rapid  dynamic  response  can  be  further  promoted  by  boiler  designs 
which  minimize  draft  losses  and  air  and  gas  side  volume. 

The  backlash  term,  a  relatively  hi^  frequency  component,  is 
caused  by  the  mechanical  lost  motion  in  the  blower  throttle  valve 
linkage. 

b.  As  discudsed  in  the  Results  section,  the  air  flow /blower 
speed  transfer  function  can  be  considered  a  proportional  conversion 
if  the  attenuation,  and  phase  shift  resulting  from  compressibility  effects 
are  assigned  to  the  blower  transfer  function.  In  order  to  complete  a 
closed  loop  description  Of  the  system  units  of  the  transfer  function 
must  be  expressed  in  terms  of  degrees  rotation  of  the  throttle  linkage, 
and  the  appropriate  dimensions  associated  with  actuator  displacement 
included.  This  breakdown  is  shown  in  the  block  diagram  of  Figure  17. 

4.  Transfer  functions  for  the  main  boiler  energy  and  mass  re¬ 
lationships  are  discussed  in  the  same  order  as  presented  in  the  previous 
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section.  Results. 

a.  The  steam  pressure /heat  input  transfer  function  for 
the  boiler  can  be  applied  to  steam  pressure  at  either  the  superheater 
outlet  or  at  the  steam  drum,  since  the  pressure  changes  are  small  and  •  - 
the  steam  flow  was  constant.  According  to  heat  balance  considerations 
the  Bode  diagram  should  indicate  a  response  characteristic  of  pure 
integration.  As  illustrated  in  Figure  11,  however,  the  transfer  function 
was  found  to  consist  principally  of  a  first  order  lag  with  low  natural 
frequency.  The  discrepancy  is  resolved  with  the  recognition  that  the 
enthalpy  of  saturated  steam  is  not  constant  with  the  saturation  pressure. 
An  increase  in  saturation  pressure  required  additional  heat  absorbtion 
to  maintain  equilibrium.  This  constitutes  a  form  of  negative  feedback 
in  the  integration  loop,  which  converts  the  integration  to  a  first  order 
lag. 

A  further  consideration  is  the  effect  of  steam  pressure  changes 
in  the  steam  output.  During  test  runs  the  steam  valve  and  all  other 
restrictions  in  the  steam  piping  constitute  a  fixed  system  resistance. 
Pressure  changes  therefore,  produce  steam  flow  changes  in  accordance 
with  critical  flow  considerations  wherein  steam  flow  varies  directly 
with  absolute  steam  pressure. 

By  the  same  argument,  va;;iation8. in  heat  input  produce  corre¬ 
sponding  variations  in  steam  temperature  with. the  result  an  added 
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variation  in  heat  output. 

Attempte  to  modulate  eteam  flow  to  maintain  constant  heat 
output  are  considered  highly  impractical.  The  transfer  function 
therefore  includes  all  three  of  the  aforementioned  effects. 

b.  The  transient  response  of  steam  pressure  to  heat 
output  is  shown  in  Figure  12.  This  response  appears  to  contain  two 
parallel  connected  transfer  functions,  a  first  order  lag,  and  an  in¬ 
tegration.  In  the  light  of  the  discussion  in  (a)  above,  it  seems  diffi¬ 
cult  to  justify  the  existence  of  the  integration  term.  Apparently  the 
transient  response  was  not  carried  to  ultimate  completion,  with  &e 
result  the  loss  of  low  frequency  data  which  might  otherwise  have  re¬ 
vealed  the  presence  of  a  linear  lag  in  lieu  of  a  pure  integration.  On 
the  other  hand,  the  integration  is  in  accordance  withithe  theory,  the 
exponential  characteristics  being  explained  by  the  side  effects  described 
in  (a)  above.  If  we  assign  the  dimensions  steam  pressure /heat  input 
and  steam  pressure/heat  output  to  these  two  transfer  functions,  only 
the  higher  frequency  responses  should  be  considered.  In  any  case, 
these  two  functions  should -be  compatible  in  form.  With  this  in  mind, 
and  neglecting  feedback  effects,  the  steam  pressure/heat  input 
tramsfer  function  can  be  modified  to  yield  the  following: 
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P/  =  15  X  10 

^Wo 


4 


_J _ 

1  +  .4p 


) 


psi 

lbs /hr 


c.  Steam  temperature  characteristics  of  the  superheater 
were  adequately  determined  by  experimental  frequency  response 
techniques.  Unfortunately,  no  linear  analytical  transfer  function  could 
be  found  to  approximate  ..the  experimental  data.  On  the  other  hand, 
th.8  response  dees  not  constitute  an  element  oi  tVie  closed  loop  system, 
and  is  therefore  cf  less  analytical  importance.  .  The  Bode  diagram  of 
Fieu.re  13  may  be  useful  in  establishir:g  criteria  for  superheater  designs 
!,".d';..'  adverse  ciansients  or  thermal  fatigue  conditions,  for  temperature 
t'-a'..sm..;cer  or  indicator  design,  etc. 

d.  The  Bode  plot  for  response  of  steam  temperature  to 
s-eam  f  ,ow  may  he  obtained  by  shifting  ihc  amplitude  curve  of  Figure  13 
di  rtV.  <'.'.5  db,  as  discussed  under  paragraph  3  d  of  "Results". 

e.  Bciiler  drum  water  level  response  to  changes  in  fuel  rate 
are  ir.dijated  in  Figure  14.  This  response  takes  the  form  of  two  series- 

:;,ed  first  order  lags.  This  ef.fecc  is  ana’.ogous  to  the  "shrink  and 
swell"  phenomenon  commonly  observed  as  response  to  changes  in 


steam  rate. 


f.  Water  level/heat  output,  or  more  properly,  water  level/ 
steam  Ucw,  was  determined  by  transieni  response  methods.  In 
addition,  to  the  linear  lag  representing  shrink  and  swell,  the  integration 
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term  obeerved  wae  a  measure  of  the  loss  (or  gain)  of  weij^t  of  contents 
of  the  boiler  due  to  differences  between  input  and  output  flow  rates. 

In  terms  of  drum  water  level,  the  rate  of  integration  (as  determined 
from  the  transient  response  curve  of  Figure  15)  was  — 1.15.x  lO'^inches 
per  second  per  pound  per  hour  steam  flow/water  flow  differential. 

g.  The  water  level/water  flow  transfer  function  obtained 
from  the  frequency  response  results  shown  in  Figure  16  exhibited  a 
low  natural  frequency  integration  as  the  primary  component.  This 
integration  results  from  the  mass  continuity  equation,  and  in  theory 
is  the  complement  of  the  water  level/steam  flow  integration  discussed 
in  (f)  above.  The  natural  frequency  of  this  function  occurs  where  the 
gain  is  unity  (zero  db),  and  from  the  natural  frequency  the  integration 
rate,  2T(fnf  found  to  be  1. 13  x  10  ^  inches  per  second  per  pound 
per  hour  steam  flow /water  flow  differential.  Considering  the  fact  that 
the  steam  flow  and  water  flow  effects  on  boiler  water  level  were  obtained 
by  different  methods,  using  different  instrumentation,  and  in  separate 
test  runs,  the  close  agreement  between  these  complementary  functions 
is  considered  nothing  less  than  remarkable. 

h.  The  steam  pressure /water  flow  transfer  function  was 
derived  analytically  and  the  natural  frequency  based  on  other  experimental 
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results.  This  analysis  appears  in  paragraph  3(h)  of  "Results". 

Further  experimental  justification  for  the  transfer  function  obtained 
is  contained  in  the  following  section. 

5.  The  closed  loop  block  diagram  appears  in  Figure  17.  Con¬ 
sideration  of  the  elements  comprising  this  diagram  gives  rise  to 
the  following: 

a.  The  energy  balance  across  the  boiler  indicates  that 
all  integration  components  in  the  heat  balance  transfer  functions 
must  add  to  zero  at  steady  state  steaming  conditions.  Based  on  a 
steaming  rate  of  56,000  pounds  per  hour,  with  fuel  rate  4200  pounds 
per  hour,  the  natural  frequencies  of  the  transfer  functions  of  Figures  11 
and  12  and  paragraph  3(g)  of  "Results  '  can  be  converted  to  Btu' s  and 
a.gebraically  summed  as  follows; 

Ws'j  =  •  .  91  X  10  ”  4  X  56,  000  =  -5.  10 
^  ^-228  X  10"'*  X  56,000  =  -i.zg 

Wq)  =  15.0  X  10  -4  ^  ^200  =  +6.30 

The  algebraic  sum,  —.08,  represents  an  experimental  error  of 
(.08/6.38),  or  about  1.25  percent.  This  error  is  considered  within 
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the  limitatlone  of  teet  inatrumentation  and  technique.  This  level  of 
accuracy  further  juatifiee  the  analytical  transfer  function  obtained 
in  paragraph  3  (h)  of  "Reaulta". 

b.  The  maaa  balance  acroas  the  boiler  alao  requires 
that  all  Integration  components  add  to  zero  at  steady  state.  Thus 

the  natural  frequencies  of  the  water  level/steam  flow  and  water  level/ 
water  flow  transfer  functions  from  Figures  15  and  16  may  be  summed 
as  follows: 

2Tr{f„l^-1.15  X  10"^ 

2T\{£n\.^^-  1.13  X  10’^ 

The  algebraic  sum,  -.02  x  10*^,  reflects  an  experimental  error, 
of  about  1.8%,  a  remarkably  high  degree  of  experimental  accuracy. 

c.  Several  of  the  elements  comprising  the  closed  loop 
systems  shown  in  Figure  17  were  not  evaluated.  These  elements  in¬ 
clude  the  Copes -Vulcan  feedwater  control  system,  including  the 
thermo-mechanical  transmitter)  the  Bailey  KP-13  steam  pressure 
transmitter;  the  General  Regulator  ratio  relays;  and  the  General 
Regulator  low  pressure  signal  selector.  In  the  case  of  the  transmitters 
the  inputs  could  not  be  independently  controlled;  in  the  case  of  the 
mechanical  combining  relay  of  the  feedwater  regulator  the  three 
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element  inpate  coald  not  be  eimalated  in  open-loop;  the  ratio  relays 
were  considered  to  have  flat  frequency  responses  in  the  range  of 
interest,  and  were  therefore  gain  devices  only;  the  low  pressure 
signal  selector  is  essentially  a  binary  logic  device,  and  its  frequency 
response  is  not  descriptive  of  its  function  in  an  analog  loop. 

The  only  devices  listed  above  that  contain  static  non-linearities 
are  the  feedwater  regulator  components.  The  thermostat  output  as 
a  function  of  water  level  in  the  boiler  drum  is  shown  in  Figure  18; 
the  flow  elements  of  the  combining  relay  produce  inputs  to  the 
controller  that  vary  with  the  square  root  of  the  flows.  These  functions 
are  included  in  the  block  diagram,  Figure  17. 

d.  The  closed  loop  frequency  response  can  be  obtained 
from  a  Black-Nichols  plot  of  the  open  loop  frequency  responae  of  the 
system.  Such  a  diagram  appears  in  Figure  19.  This  diagram  re¬ 
presents  the  response  of  the  closed  loop  to  disturbances  in  the  set 
point  steam  pressure.  From  the  Black-Nichols  diagram  it  is  evident 
that  the  maximum  gain  of  theclosed  loop  is  9  db,  occurring  at  a 

frequency  of  .004  cycles  per  second.  The  closed  loop  phase  angle 
o 

is  -60  .  The  closed  loop  diagram  indicates  that  stability  criteria 
are  not  being  met.  Response  to  an  increase  in  set  point  indicates 
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excessive  low  frequency  phase  shift,  that  to  a  decrease  in  set  point, 
insufficient  loop  g^in...  These  two  conditions  bracket  the  ideal  closed 
loop  gain  and  phase  stability  criteria,  indicating  that  the  system  is 
optimally  adjusted  while  at  the  same  time  inherently  incapable  of 
optimum  performance.  These  conclusions  are  supported  by  observa¬ 
tion  of  control  system  performance  on  the  Laboratory  boiler.  From 
the  diagram  it  is  further  obvious  that  one  method  for  improvement 
of  the  system  response  lies  in  elimination  of  the  reset  mode  in  the 
steam  pressure  control  loop.  This  in  fact  is  exactly  the  technique 
employed  in  the  systems  for  DD945  and  OLG-9  Class. 
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